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I Abstract. The masses of charmonium states immersed in nuclear matter are calculated in LO QCD 

. and in QCD sum rules. While the mass shift for7/i// are found to be less than -10 MeV, those for the 

Xq.1.2 and i//(3686) and i//(3770) are found to be more than -40 MeV. We investigate the feasibility 
Q [ of observing such mass shifts in the future accelerator project at GSl. 

o 

^ : INTRODUCTION 
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^ ■ Understanding hadron mass changes in nuclear medium and/or at finite temperature can 

O ! provide valuable information about the QCD vacuum|[ll0,|3]. While the mass shifts for 

hadrons made of light quarks are sensitive to the restoration of the spontaneously broken 
■ chiral symmetry breaking 1 2!, those for the the heavy quark systems are sensitive 

to the changes of the non-perturbative gluon fields in nuclear matter. For the J/VL which 
consists of a charm and anticharm quark pair, both the QCD sum rules analysis f/, oj and 
the LO perturbative QCD calculation [9, 10] show that its mass is reduced slightly in the 
nuclear matter mainly due to the reduction of the gluon condensate ((^G^)) in nuclear 
matter, which is expected to decrease by 6% at normal nuclear matter density. However, 
^ I the changes are much larger for excited charmonium states, due mainly to larger color 

C ■ dipole size of these excited states. 

In this report, we summarize the expected mass shift for charmonium states in nuclear 
kJ \ matter and study the feasibility of observing such mass shift in the future accelerator 

^' project at GSlflT]. 

The lowest dimensional QCD operators that characterizes the non perturbative nature 
of the QCD vacuum are the quark and gluon condensate. These condensate are estimated 



m 

O 



to have the following large non-perturbative expectation values in the vacuum Ill2|] . 



i^F^,) ~ 1.5GeV/fm3, 

(qq) ~ 2fm-3. (1) 

The gluon condensate can be written as the difference between the magnetic = F^- 

and electric = ^F^j condensate, which respectively contribute to half of the zero 
temperature gluon condensate, 

{^B^) = -{^E') = U^Fi,). (2) 
n K z n ^ 



The above relation follows naturally from the Euclidean formulation of QCD at zero 
temperature, such as in the lattice QCD, where the Euclidean space electric (mag- 
netic) condensate is defined with a minus (plus) sign relative to its Minkowski space 
counterpart[12]. Due to the symmetry in the time and space directions in the 4 di- 
mensional Euclidean space, the Euclidean space electric condensate is expected to have 
the same expectation value as the magnetic one l.l3Lil4.] . Hence the relation among the 
Minkowski space condensate in Eq.© follows. 

At nuclear matter, the non perturbative quark and gluon field configuration are ex- 
pected to change appreciably, such that the average gluon and quark condensate values 
decrease by 6% and 30%, respectively. These model independent results are obtained 
from the linear density approximation and the nucleon expectation values of the quark 
and gluon condensate, which are respectively known from the experimentally measured 
pi-N sigma term and from taking the nucleon expectation value of the trace anomaly 
relation! 15]. The electric and magnetic part of the gluon condensate at nuclear matter 
can be estimated separately, by using the twist-2 gluon operator, 

{N{p)\y2rF^Fav\N{p)) = {p^py - ^m^g^y^ 2A2(g), (3) 

where A2(^) is the second moment of the gluon distribution in the nucleon and is around 
0.45, when the renormalization scale is in the of order 1 to 2 GeV. Using this and the 
linear density approximation, we have, 

2v /4 0,3 2 . \Pn.m. 

= -[^^M^M-^nr,-A2j^^. (4) 

Here, ~ 0.75 GeV is the mass of the nucleon in the chiral limitlH^, which comes 
from taking the nucleon expectation value of the trace anomaly relation = — f^^^v 
As can be seen from Eq.©, due to the additional factor of ^ in the second terms, the 
changes are dominated by the contribution from the first terms. 



CHARMONIUM MASS SHIFT FROM QCD 

The mass shift of charmonium states in nuclear medium can be evaluated in the pertur- 
bative QCD when the charm quark mass is large, i.e., nic — 0°. In this limit, one can 
perform a systematic operator product expansion (OPE) of the charm quark-antiquark 
current-current correlation function between the heavy bound states by taking the sep- 
aration scale (jU) to be the binding energy of the charmonium llSl l. The forward 
scattering matrix element of the charm quark bound state with a nucleon then has the 
following form: 

T{q^ = m\)=Y^-^^{ff,)^. (5) 



Here, Cn is the Wilson coefficient evaluated with the charm quark bound state wave 
function and {^n)N is the nucleon expectation value of local operators of dimension n. 

For heavy quark systems, there are only two independent lowest dimension operators; 
the gluon condensate ((^G^)) and the condensate of twist-2 gluon operator multiplied 
by Us ({^GajiGy)). These operators can be rewritten in terms of the color electric and 
magnetic fields: {^E^) and {^B^). Since the Wilson coefficient for {^B^) vanishes in 
the non-relativistic limit, the only contribution is thus proportional to {^E^), similar to 
the usual second-order Stark effect. We shall thus calculate the mass shift of charmonium 
states due to change of the gluon condensate in nuclear medium by the QCD second- 
order Stark effect 113]. 

The mass shift of charmonium states to leading order in density is obtained by 
multiplying the leading term in Eq.©, by the nuclear density p]\/. This gives. 



Am^(e) = / dk^ 



^eA (6) 

In the above, rrif^ and are the nucleon mass and the nuclear density, respectively; 
{^E^)i\i ~ 0.5 GeV^ is the nucleon expectation value of the color electric field obtained 
from Eq.© and £ = Inic — ruxf,. In Ref.L9], the LO mass shift formula was derived in 
the large charm quark mass limit. As a result, the wave function ^f{k) is Coulombic 
and the mass shift is expressed in terms of the Bohr radius uq and the binding energy 
£o = Im-c — nij/yf. This might be a good approximation for 7/i// but is not realistic for 
the excited charmonium states as Eq.® involves the derivative of the wave function, 
which measures the dipole size of the system. We have thus rewritten in the above the 
LO formula for charmonium mass shift in terms of the QCD parameters = 0.84 and 
Mc — 1.95, which are fixed by the energy splitting between J/xjr and i//(3686) in free 
spacef^. Furthermore, we take wave functions of the charmonium state to be Gaussian 
with the oscillator constant /3 determined by their squared radii (r^) = 0.47^, 0.74^, 
0.96^, and 1 fm^ for Jjjr, Xo,i,2, V^(3686), and v/(3770), respectively, as obtained from 
the potential models [19]. This gives /3 = 0.52, 0.43, 0.39, and 0.37 GeV if we assume 
that these charmonium states are in the 15, IP, 25, and ID states, respectively. Using 
these parameters, we find that the mass shifts at normal nuclear matter density obtained 
from the LO QCD formula Eq.® are -8, -40, -100, and -140 MeV for 7/1//, XoX2' 
V/(3686), and V/(3770), respectively [20]. 

Although the higher twist effects on the charmonium masses are expected to be 
nontrivial, the result for 7/i// is consistent with those from other non-perturbative QCD 



studies, such as the QCD sum rules [|7|,|8|] and the effective potential model 112 1LI22L 12311 . 
which are all based on the dipole interactions between quarks in the charmonium and 
those in the nuclear matter. The QCD sum rule results can also be applied for the Xo.i.2 
states, and the results from the leading order gluon condensate is summarized in Table. 

ID 

Higher twist effects can be estimated in some calculations. For QCD sum rules for 
7/1//, the corrections coming from dimension 6 operators are less than 30% of the 



leading order resultsfs']. The contributions from the DD meson loops in the 1^(3686) 
and v/(3770) are also found to be less than 30% of the LO QCD resu1tL2nll. 
All the results are summarized in table[T] 



TABLE 1. Charmonium Mass shift in nuclear matter in MeV 



Charmonium 




QCD 2nd order Stark Effect 


QCD sum rules 


Effects of DD loop 




0-+ 


-8 MeV 


-5 MeV 


No effect 




1— 


-8 MeV 


-7 MeV 


<2MeV 


Z0,l,2 


0,1,2++ 


- 40 MeV 


-60 MeV 


No effect on Xi 


Va(3686) 


1 — 


-100 MeV 




< 30 MeV 


1/^(3770) 


1 — 


-140 MeV 




< 40 MeV 



OBSERVABILITY 



Since the mass shift of the heavy quark system reflects the changes of the Gluon 
field configuration in the vacuum, it would be interesting to observe such effects in 
experiment. 

Consider an anti-proton with incoming four momentum (ft), 0, k) annihilating a proton 
at rest (mjv,0,0) and creating a charmonium moving with velocity v. The required 
incoming momentum k to create a charmonium state are summarized in Table I. 



TABLE 2. Required momentum to create Charmonium with outgoing 
velocity v 







J/W 


Xo 


Xi 


Xi 


V/(3686) 


V/(3770) 


k (GeV/c) 


3.7 


4.1 


5.2 


5.5 


5.7 


6.2 


6.5 


(0 (GeV) 


3.8 


4.2 


5.3 


5.6 


5.8 


6.3 


6.6 


v(c) 


0.78 


0.8 


0.83 


0.84 


0.85 


0.86 


0.87 



Hence, the required incoming energy of the anti-proton to produce the charmonium 
state range from 4 to 6 GeV. In these energy region, the absorption cross section 
Op^p ~ 50 mb. Hence, the anti-proton would be absorbed after travelling less than 1 fm 
in the nuclear matter. Moreover, once a charmonium state is created, their speed would 
be less than 0.9 c, which means that it will have to travel more than 10 fm/c to pass the 
diameter of a nucleus of A=125. Hence, considering the increased width of charmonium 
due to nuclear absorption0|, the charmonium is expected to decay inside the nucleus. 

The cross section for the production of charmonium states and its subsequent decay 
into dileptons or 7/i// + 7 states are given by the following Breit-Wigner formula 



where k is the cm. momentum, E is the cm. energy. Bin and Bout are the branching 
fractions of the resonance into the entrance and exit channels. The 2s + I are the 
spin multiplicities of the incident spin states and J the spin of the charmonium. Also 

^ Total —medium ^ Total "t~ ^medium • 



If we substitute the medium mass shift and increase in width (due mainly to collision 
broadening), the cross sections are in the order of one to few hundred pbarn 

The expected luminosity at the anti proton project at GSI is 2 x lO^^cm^^s^^ . There- 
fore if the cross section is Ipb, it would corresponds to about 17 events per day. 



TABLE 3. Measurable decay channel and expected event rate at GSI future accelerator. 



Charmonium 


^Total ~l~ ^medium 


Final state 


cross-section to final state 


events per day 




y/VA(3097) 
1/^(3686) 
V/(3770) 


87 KeV+ 20 MeV 
300 KeV+ 20 MeV 
23.6 MeV +20 MeV 


e+ + 


6pb 
0.6 pb 
1 pb 


100 
10 
17 


Zco(3417) 
Zci(3510) 
Zc-2(3556) 


16.2 MeV+ 20 MeV 
0.92 MeV+20 MeV 
2.08MeV+20 MeV 


J/w+y 
J/w+y 
J/w+y 


200 pb 
80 pb 
350 pb 


3400 
1360 
5950 



Hence, the mass shift will be observable in the anti-proton project at the future 
accelerator facility at GSI. 
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